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Abstract 

The variety of the phase transitions in Induced QCD are studied. Depending upon 
the parameters in the scalar field potential, there could be infinite number of fixed 
points, with different critical behavior. The integral equation for the density of the 
eigenvalues of the scalar field are generalized to the weak coupling phases, with the 
gap at the origin. We find a wide class of the massive solutions of these integral 
equations in the strong coupling phases, and derive an explicit eigenvalue equation for 
the scalar branch of the mass spectrum. 
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1 Introduction 



Induced QCD was suggested as a possible model of hadrons two months ago IJ. This 
is the lattice model of the scalar field $ in adjoint representation of the SUn gauge group, 
interacting with the usual lattice gauge field Ui, defined at links I. The unusual feature, which 
allows one to solve the model exactly in the large N limit [J for arbitrary dimension D 
of the lattice, is the absence of the bare self-interaction for the gauge field. 

The idea is, that this self-interaction would be induced at larger distances, where the 
scalar field decouples, being heavy. Such induction often takes place in two-dimensional 
models of gauge and gravitational fields, such as CPjv models, or the Liouville theory. 

The parameters of the scalar potential U(*f>) should be carefully adjusted for this miracle 
to occur in more than two dimensions. The bare mass = U"(0) should be chosen so, that 
the effective scalar mass becomes much less than the lattice cutoff. This effective mass 
m e ff would serve as the ultraviolet cutoff for perturbative QCD. 

At the scales less than this mass, the scalar field could be integrated out, which yields the 
Yang-Mills term N/3 trF^ v in effective Action, along with a variety of unwanted higher order 
terms. The bare quartic coupling should tend to a hypothetical critical point, to provide a 
large positive coefficient (3, which could serve as the bare coupling of perturbative QCD. Then 
the physical mass scale m phys would automatically come out m phys ~ m eff exp (— A/3) /3 B , 
with coefficients A, B known from perturbative QCD. 

The computation of induced gauge coupling (3 as a function of the bare scalar mass and 
the bare quartic coupling A = U""(0) is not an easy task. Apparently, the quartic coupling 
should tend to some ultraviolet fixed point of the renormalization group, corresponding to 
both gauge and scalar charges. The recent computer simulations M for the SU2 model 
showed the line of the first order phase transitions in the m.g, Ao plane, with the endpoint, 
which might correspond to the second order phase transition. 

The general equations of were studied in great detail for for special cases D = 1 
and quadratic potential U{4>) ||. In agreement with computer simulations and theoretical 
expectations, there is no critical point of the QCD type for the quadratic potential. At 
D = 1 in continuum limit the familiar fermionic solution was reproduced. The more general 
solution, for the finite lattice spacing at D — 1 was found in |J , again in complete agreement 
with previous work on matrix models. 

There are some lattice artifacts in this model, which it is not yet clear how to remove. 
Namely, there is an extra local symmetry, which must break to allow quarks to move 
inside hadrons. The general mechanisms of this symmetry breaking were found long ago by 
Khokhlachev and Makeenko J7|, and later confirmed by numerical experiments and mean 
field analysis [|, 

It was recently conjectured in [12], that the same mechanisms spontaneously break Zjv 



symmetry in this model at the critical point, which corresponds to the continuum limit of 



the lattice theory. In a subsequent paper |T3 it was argued on the basis of the mean field 



analysis, that the transition must take place before the critical point, i.e. still in the 
strong coupling phase of our model, if it really induces QCD. 



Various properties of the generalized Wilson loops in this model were studied in [14 
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Some interesting mathematical structures were found, which could be used in gauge and 
string models regardless conjectured induction of QCD. 

There is the general heuristic argument , that the nontrivial fixed point of this model 
could be nothing but QCD. The argument is based on the common belief, that there is only 
one nontrivial theory in four dimensions: the asymptotically free, quark confining QCD. 
Even if this argument fails, the confining solution of Induced QCD would be an exciting 
alternative to the usual QCD. This would be the first solvable model of the QFT in four 
dimensions. 

In the second paper we came very close to this goal, by reducing the solution of the 
model to the following nonlinear integral equation for the vacuum density p(A) of the scalar 
field eigenvalues 



RW = ^u'(\) + ^p I dA'f^, (1.2) 



p 

where 

f ' lX) = 2ZT vv ' D ^J-oc ""' A - A' 

and U(4>) is the scalar field potential. We call this equation the master field equation, or 
MFE, because this density could be regarded as the long-sought master field of QCD. We 
found exact powerlike solution at A — > 00 

p(A) ex |A|°; cos?ra = ~ 3J ^_ g i a > 1, (1.3) 

which was quite encouraging, since the scaling index a showed no pathologies, like those of 
the string models. 

The forbidden interval here is | < D < 1, and the only rational values are a = n + \ 
at D = 0, a = 2n at D = ~, a = 2n + 1 at D = 1 and a = 2n + 1 ± | at D = 2. The 
solutions at D = 0, | are unphysical, as the assumption of vanishing density at the origin is 
never satisfied in the corresponding matrix models. D = corresponds to the one matrix 
model, and D = ~ correspond to the two matrix model (the number 2D of links meeting in 
each cite equals 1 here, as there is only one link, connecting two cites). D = 1 solution is 
already physical, as we discuss in more detail below. The |A| 2n+1 singularity comes about 
as the singularity of the tip of the upside-down even potential in the usual solution of the 
D = 1 matrix models. As for the first nontrivial case D = 2, unfortunately, the adjoint 
scalar field model cannot be solved by conventional methods even at D = 2, so there is 
nothing to compare with this solution. 

The jj expansion was considered in the third paper ||, where we found the integral 
equation for the propagator of the effective field theory with p(A, x) as dynamical field. This 
linear equation involves the vacuum density in its kernel, and for the powerlike density, the 
powerlike solutions for the corresponding wave functions in A space were found. 



3 The eigenvalues A have dimension m^ D 1 so that at D > 2 the physical region is A — > in the lattice 
units we are using. 



2 



Still, the solution is incomplete, as there is no mass scale. This is the solution exactly 
at the critical point, where there are scaling laws in the A space. The physical solution of 
Induced QCD must involve the mass scale, which requires the more general solution of the 
MFE. 

In this paper we find an infinite family of such solutions, which turns out to be a particular 
superposition of the previous powerlike terms. The implications of this simple observation 
are very interesting. Now, there is a calculable mass spectrum with nontrivial scaling indices 
in arbitrary dimension D > 1. 

We generalize the Riemann-Hilbert problem for the weak coupling phase, where there is 
a gap at the origin in the density of eigenvalues. We present the new derivation, which, as 
we hope, is easier to comprehend, than that of 0. 



2 Massive Solution of the Riemann-Hilbert Problem 

The classical equation ( |i.l| ) in the local limit, when r(A) = R(X) — X ~ p(A) <C A was reduced 
in the previous paper to the following nonlinear boundary problem. 
Let us introduce two functions 



/DO P 2 (/^) 
dp, . (2.2 
-oo /i — Z 



At z — > X + tO , we have 



V(z)^ T ^r(X)+mp(X). (2.3) 

Here U((p) is the bare potential, but the contributions from the large eigenvalues ( of the 
order of the lattice cutoff) in the integral effectively renormalize this potential. In the local 
limit, at A <C 1 in lattice units, the density p ~ |A| a , so that formally the integral diverges 
if we substitute the local density, i.e., it is dominated by the lattice scales A ~ 1, where 
the solution is not universal. The corresponding number of subtractions should be made, 
as usual in dispersion relations, or, which is more convenient, one could make the analytic 
continuation of these integrals in the scaling dimension a from the convergence region a < \ . 
The subtraction polynomial renormalize the bare potential U. 

Both functions are analytic in the upper half plane and have the symmetry property 

V(-z) = -Viz); Q(-z) = -Q{z). (2.4) 

In other words, real (imaginary) parts are odd(even) with respect to the real part of z. 

As was discussed already in the first paper, the renormalization group analysis tells us, that the log- 
arithmic laws of the asymptotic freedom translate into the power laws in the induced QCD models. The 
critical indices depend upon the effective quartic scalar interaction, which is not calculable by perturbative 
methods. 
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At ^z — > +0 by construction 

SQ = (^p) 2 . (2.5) 

On the other hand, as shown in in virtue of the classical equation up to 0(V 3 ) 

termsf] 

m = (p 2 ) . (2.6) 

The last three equations represent the nonlinear Riemann-Hilbert problem. 
In the previous papers only the massless solutions were found 

V = ^^C, (2-7) 
cos ^ 



,a(-^) 



2d 



Q = '—. (2.8) 

cos7ra 

Various values of a = In + 1 ± ^ arccos 3J ^_ 2 correspond to various fixed points of the model. 

We could not find any general theory of the nonlinear Riemann-Hilbert problem, but in 
this particular one we found the class of exact solutions, which are built from the above 
power terms, 

cos7ra 

It is not difficult to check this solution. The symmetry property is manifest, and so is 
the first equation. As for the second equation, the key point is that at z = A + iO 

= A\X\ a + B\Xf; m> = ^tan^ (A\X\ a - B\\f) , (2.11) 

|A| 2 

so that the cross terms in ^s(V 2 ) = 23'PK'P are absent. 

In general, these A and B could be arbitrary polynomials of z 2 . This ambiguity reflects 
the ambiguity in the choice of the initial scalar potential U((p). However, the critical phe- 
nomena, which we are interested in, are universal, as they take place in the infinitesimal 
vicinity of the origin, in the lattice units we are using. 

In this limit, only the two leading terms in A, B can be left. With proper redefinition 
of a, (3 this corresponds to constant A, B. We assume, that a > (3, then the critical region 
corresponds to 

,~» S (|)" ;w =-i^, (2.12) 

so, that B = at the critical point. The leading index a must be greater than 1, as it follows 
from original derivation. As for the subleading index /3, it should be greater than —1 for the 
density to be integrable at the origin. 



These terms are down by a power of the ultraviolet cutoff in the local limit. 
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Using the language of renormalization group, the /3-term represents the perturbation of 
the UV-stable fixed point a by the relevant operator of the lower scaling dimension. In 
general, the coefficient B linearly vanishes at the critical point. At B < the solution 
becomes unstable, as the density changes sign near the origin. In this case, there would be 
the phase transition to the weak coupling phase. The local limit of the MFE is different in 
this phase, as we shall see in the next Section. 

The lowest solution for a, (3 in the strong coupling phase would be 1 < a < 1.5, and 
(3 = 2 — a. In four dimensions 

{a, (3} = {1.36901,0.63099}. 

The next one is 

{a, (3} = {2.63099, -0.63099}. 

Note that the singularity in p is still integrable for this solution. These 
solutions for k — 1. At the next level, k = 2, there are three solutions 

{a, (3} = {{2.63099, 1.36901}, {3.36901, 0.63099}, {4.63099, -0.63099}}. (2.15) 

We discuss the choice of the solution later, when we study the wave equation. As we shall 
see, the first solution at k — 1 have tachyons. 

It is worth mentioning, that at D = 1 the equations degenerate. The real part of Q 
vanishes, so that the most general solution with proper symmetry would be iW(z 2 ), where 
W is some polynomial with real coefficients. Then, from the first equation we find 

np(\) = ^W{\ 2 ). (2.16) 

This is in complete agreement with the solution of the lattice MFE, found recently by 
D. Gross ||. In this case W(z 2 ) = 2tx 2 {E — U(z)) where U(z) is initial potential, and the 
chemical potential E is to be determined from the normalization of density. This solution 
also agrees with conventional solution of the D = 1 model in terms of effective fermi gas. 
Note, that the generic singularity is |A| 2n+1 , for W(X) oc A 4n+2 . This agrees with above 
powerlike solution. 

With our "S-matrix" approach we immediately find the solution, but cannot relate the 
parameters to those of original lattice theory. Fortunately, this is never needed. What is 
really needed, is to check internal consistency of the solution, such as absence of ghosts and 
tachyons, which is not apriori guaranteed in the "S-matrix" approach. 

Another comment. At any D there always exists a trivial solution, without critical 
behavior 

V = za{z 2 ) + i^ 2 ); Q = %b 2 {z 2 ) + 2^-^za{z 2 )^fb(^), (2.17) 

with the support of eigenvalues at b(\ 2 ) > 0. In the simplest case of constant a and linear 
b this is the semicircle solution. As was recently shown in pL for the case of quadratic 
potential U this is the only solution in the strong coupling phase. 

Unfortunately, the simplest nontrivial scaling solution would take at least two adjustable 
parameters, so the higher order terms in a potential are required. In this case, as we suspect, 
the explicit solution of the lattice model is unavailable, and one either has to rely upon the 
above Riemann-Hilbert approach, or use the numerical methods to solve the lattice MFE. 



(2.13) 
(2.14) 

are the only physical 
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3 Phase Transition 



The Riemann-Hilbert problem was derived under assumption, that there was infinite support 
of the eigenvalues, without any gap at the origin. As was mentioned in the first paper [|TJ, 
we expect this model to undergo the phase transition from above strong coupling phase to 
the weak coupling phase, with the gap in the support of the eigenvalues. In that paper we 
could not find nontrivial spectrum, because the kinetic term in the effective Lagrangean for 
the density fluctuations vanished at N = oo. 

The more recent solution || produces such term regardless the phase of the model. 
The term proves to be positive definite, which means that this solution is different from 
the first one. Perhaps, there was something wrong with the assumptions of the orthogonal 
polynomial method in this case.f] 

Anyway, let us assume, that there is a gap from —a to a in the vacuum density p(A). 
The basic equation (|1 . 1| ) remains the same, but the dispersion relation between real and 
imaginary parts of all analytic functions modifies. The simplest way to account for these 
changes is to note that conformal transformation 



C(z) = Vz 2 - a 2 (3.1) 

maps the upper half of z plane onto the upper half plane of C, eliminating the gap from —a 
to a. 

We could use old dispersion relations for the even functions of z with C instead of z. For 
the odd functions there would be fictitious singularity at z 2 = 0. In particular, the odd 
function Viz) was reconstructed in from the real part. The derivative of this relation 
reads (so far, at a = 0), 

' [Z) 2(1 -D) + J-oo vr \z - C(u) z-R(u))' [ ' 

where 

C(v) = R(v)-vk P {v). (3.3) 

This is the difference of two Cauchy integrals, the first one going over the complex curve 
z = C(v) in the lower half plane, and the second one going backwards over the real axis 
z = R{y). One could write this as a single Cauchy integral over the complex contour 
C = {C,R}, 

w/ . U"(z)-2D r dyWAy) 

V (z) = — H v- +i <P — — , (3-4) 

K J 2(1 - D) JC 7T z - y ' K J 

where the density Wq is given by the parametric equation 

Wc (C(u)) = ^- (3.5) 



6 In some cases, the even and odd coefficients of expansion in orthogonal polynomials tend to different 
analytic functions of jr. Certainly, this issue must be further analyzed. 
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Transforming this relation from z to £, we find 



w 2(1 - D) Jc 7T CO) -y 

Wc(C(C{u))) = -^ (3.7) 

In terms of original variables 

V(z) = U " {Z) ~ 2D + i t~ - ( aC{v)) - Cim) ) (3 8) 

1 ] ~ 2(1 -D) + loo 7r UW-C(C(»/)) CO)-C(^))/ 1 ' J 

Note, that the integrand vanishes inside the gap, as C{v) = R(v) there. Also, note, that 
in virtue of the symmetry R(—v) = —R(u),p(—u) = p(u) the real (imaginary) part of 
this function is odd(even) function of 3fe, as it should be. One can readily check that this 
function is analytic in the upper half plane, that its imaginary part at the real axis vanishes 
at z 2 < a 2 , and that its real part at z 2 > a 2 agrees with derivative of MFE. [] 

Let us now go to the local limit in above complex contour integral for V(z), using the 
same expansion, as in ref. @, 

R(u) = v + r(i/); G{v) = u + c{u); c{y) = r(y) - inp(y) < v. (3.9) 

Iterating implicit equation for Wc, we find 

Wc(((z)) = 1 - c'(z) + \ (c 2 (z))" + 0(c 3 ) (3.10) 

We are interested in 7T//(A) = QV'(X + lO). One contribution to this imaginary part 
comes from QW(£(A)) times the ^-function term at y — C(A). Subtracting the same terms 
with p = 0, to account for the integral with C(u) replaced by R(v) , we find 

9fWb(C(A)) - 9W fl (C(A)) = np'(X) - TV (r(A)p(A))" (3.11) 

We see, that the left side exactly cancels with the first term, so that we have to keep the 
0(c 2 ) terms. 

There is the second contribution to ^sV'(X + zO), coming from the principal value of the 
integral of $lWc- As before, subtracting the integrals with C and R , we find 



-i p L* (VA--^-») W^' (3 - 12) 



or, symmetrizing in y — > —y and transforming ydy = udu, 



-TV 



^r vX^^i W- « 3 - i3) 



TT 



7 The sceptical reader is invited to check this formula by means of usual dispersion relations in the z 2 
plane, with extra factors y/a 2 — z 2 to convert real part to imaginary. In differentiating the MFE, one should 
take into account the 6- function terms, coming from the discontinuity of the actangent at ±oo. 
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where S = {(— oo, —a), (a, oo)} is the support of eigenvalues. 
Collecting the terms, we arrive at the following equation 



2(r(\)p(X))" = - X ll-^p I . w (3.14) 



2 



dv (p 2 (u))" 



X 2 ° Js.f^4 A 



At a = this equation reduces to the old one, after two integrations by parts. We drop the 
divergent polynomial terms, keeping in mind the corresponding number of subtractions in 
dispersion relations. 

The dispersion relation for r(A) reads, as before 

U'(X)-2DX D-l r , p{v) , olc , 



2D D Js X — v 

Let us now reduce these equations to the Riemann-Hilbert problem. The first function, 
V(z) is the same as before. The second function, Q(z) is introduced as follows 

e-W = x / -JtLf M (3.i6) 

This function has the same symmetry and analyticity properties, as the first one, includ- 
ing the gap in imaginary part at z = A + zO 



1 _ |! = ^ (A 2 - a 2 ) [(3P) 2 ]" (3.17) 



A 

As for the real parts, they are related at A 2 > a 2 as follows 



^l_^ Q " = i_^[c^ 2 ]" (3.18) 
At finite gap, we cannot eliminate the derivatives, because we cannot include the factor 



1 — ^2 in Q"(z) without introducing the singularity at z — 0. Still, the equations are so 
simple and universal, that one may hope to find the analytic solution, like the one we found 
in the strong coupling phase. 



4 The Mass Spectrum in the Strong Coupling Phase 



Let us now substitute the above general solution (|2.9|) into the wave equation, found in the 
previous paper JJ. In present notations, with 

~ . , DU"(z) - 2D 2 - 2D + 2 , 2 
M 2 ff (z) = D _ + polynomial = r + r x z 2 + . . . , (4.1) 
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the wave equation reads (at z = A + i0) 



(4.2) 



where P is the Euclidean 4-momentum of the vacuum excitation, corresponding to plane 
wave fluctuations of p, 

c / \ \ 1 lPx f°° , $sJ r (v + lO) 1 fA 
Sp(X,x) = -e iPx p dv \ ' -. (4.3) 

N J-oo 7T a p(U) V — X 

The pq- term in the ratio of real and imaginary parts of V yields the 5(X) term which 
drops provided 

%V = o at z = o. (4.4) 

This boundary condition selects the physical solutions. 

In the simplest nontrivial case of the mass term plus quartic interaction U (A) = |m,Q + 
|AoA 4 there are two terms To, T\ present in effective mass term. The first term To must vanish 
as Zq -1 to be relevant in the critical region z ~ Zq. This yields the equation 

r oc B s «; 5 = ^1. (4.5) 

a — [3 

In terms of the original parameters of the scalar potential, above relation describes the curve 
of the first order phase transitions in the m^Xo plane, ending at the critical point, where 
T = B = 0. This is is qualitative agreement with the simulations of [||. 

When the a > 3 solution is taken, the first t\ correction to M|^j is relevant. The similar 
estimate yields the scaling relation 



Tl oc B \ 5 1 = (4.6) 
a — p 

This would correspond to the tricritical point. In general, for a > 2m + 1 the T m z 2m terms 
in M 2 jj, coming from the <^ 2m + 2 terms in U(<p), become relevant, 

c a - 2m - 1 

r m oc B Sm ; S m = — - . (4.7) 

a — p 



Let us denote 
Note that the mass indexes 



a = k + p; (3 = k- p; p > 0. (4.8) 

. Ik — 2m — 1 

S m = - + , 4.9 

2 p 

are trivial only for the simplest fixed point, with k = l,m = 0. In general, these are 
transcendental numbers, which agrees with the induced QCD scenario, and contradicts the 
Gaussian fixed point for the scalar field. 
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Consider the infinite sum of power terms for T 

Z — t ■ wis — e) ' 

e sin \ ' 

where s = {0, 1} is the "A-parity" 

F{-z) = Sp(-X,x) = (-l) s 8p(X,x). 



(4.10) 



(4.11) 



Differentiating real and imaginary parts of V, J 7 , multiplying by and collecting the power 
terms, we find the following equation 



E / (c) A e [--(P 2 + M e 2 7/ (A)) (A\» + BA-") 
= E / ( e ) A£+fc_1 U 2 $2(e)A 2M + 2AB$ (e) + J B 2 $_ 2 (e)A~ 2 ^ 



(4.12) 



where 



Lj(fc+,i-i)j 



M e//(A) = E ^A 

m=0 

/ ?r(e- s) 
V 2 
7r(e — s) ^ 



2m 



(4.13) 



$ ±2 (e) = ecot 



± 



'(D-l) 



D t ■ (k I //)/; tan () 



$o(e) = ecot 



D 2 + 2(D-1) 



D 



tan ■ 



net 



Let us consider the simplest case k — 1, where <5 = |, and M 2 ^j = r . In this case, it is 
clear from the above equation, that /(e) = unless 



e = e — n/x; n = 0, 1, . . . , 
where the highest power e is determined by the equation 

$ 2 (e ) = 0. 



(4.14) 



(4.15) 



This highest power term was already found in the previous paper 

We solved this equation numerically in four dimensions and we found the following values 
of e for two lowest values of a A , in spectroscopic notations A = (— l) s 



1.36901+ 
1.36901" 
2.63099+ 
2.63099" 



e = {2.08496, 4.11512, 6.13072, 8.14028, . . .}, 
e = {1.05590, 3.10290, 5.12399, 7.13601, . . .}, 
e = {1.94572, 3.91608, 5.89759, 7.88500, . . .}, 
e = {2.92897, 4.90587, 8.88011, 10.87223, . . .}. 



(4.16) 



Let us pick up a particular e , and let us study the arising recurrent equation for the 
coefficients 



f{e) = y ZW k (e)f(e + kn), 



(4.17) 



k=i 
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where 



Wi(6) = -^#, (4.18) 



2A$ 2 (e) 

2B$ (e) 

' A$ 2 (e) ' 
£(P 2 + Tq) 
' 2A 2 $ 2 (e) 
£? 2 $_ 2 (e) 
' A 2 $ 2 (e) ' 



W 2 (e) = 
W 3 (e) = 
W 4 (e) = 

We could write down the formal solution 

/(6b - nil) = n ( E ^( £ o - M exp ^ /(eo ), (4.19) 

l=n \fc=l \ / / 

where it is implied, that /(e) = at e > e , and the operator ordering is as indicated, i.e., 
from I — n to I — 1. 

These coefficients should terminate at the smallest e > — (3, according to our boundary 
condition (|4.4fl 

/(e - n 0J u) = 0; n = f-^] . (4.20) 

A* 

This provides us with the spectral equation, which is a polynomial in P 2 . The roots £ n of 
this polynomial correspond to the particle spectrum. Restoring quantum numbers, 



-P 2 = T + VAB£n(s,e ). (4.21) 

The corresponding roots for the lowest levels of eo are 

£ n = ±{6.94993 ±25.62917^, 11.64012 ± 7.61785?}, (4.22) 
£ n = {0,0,0,0,0}, 
£n = 0, 

£n = {0,0,0}, 

£ n = ±49.97203, 

£n = {0,0,0}. 

With large enough r there exist tachyon-free solutions for the second fixed point, a = 
2.63099, but apparently, there are no solutions with infinitely rising masses, because of the 
sign degeneracy £ — > — £ in both fixed points. 

This is another indication of the triviality of the k = 1 fixed points (they are the only 
ones with rational critical index 5q = | of the mass spectrum). Most likely, there is the finite 
number of the free scalar particles, i.e., these are just the Gaussian fixed points. 

The case k > 1 is much more complicated, as there are also negative integer powers of A 
involved, apart from powers of A _/i in the expansion. We are going to study this case in the 
next paper. 



1.36901+ 


eo 


= 2.08496 


1.36901" 


eo 


= 1.05590 


2.63099+ 


eo 


= 1.94572 


2.63099+ 


eo 


= 3.91608 


2.63099" 


eo 


= 2.92897 


2.63099" 


eo 


= 4.90587 
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To summarize, there are two branches of the spectrum. The odd A-parity states represent 
the scalar particles, dressed by interaction with the gauge fields. The even A-parity states 
represent the mixture of "glueballs" with the even number of scalar particles. The mass 
scale behaves as certain irrational power of initial parameters of the lattice model. 

In the strong coupling phase for the two simplest fixed points with rational scaling indexes 
for masses, we computed the particle spectrum. One fixed point, with p(0) = 0, turned out 
unstable, and in the other one, with p(0) = oo, there were stable solutions. However, the 
spectrum did not rise to infinity, as one would expect in QCD, which indicates triviality of 
these fixed points. We leave for future study the exciting numerical problem of computing 
masses from above equations of higher critical points in the strong coupling phase. 

So far, we cannot even tell, whether the spectrum terminates, and whether there are 
tachyons. Perhaps, some general inequalities can be derived to answer this question. One 
way to guarantee the absence of tachyons in the given fixed point is is to arrive at this fixed 
point from the lattice MFE, with real potential, stable at infinity. This requires the serious 
numerical study of the MFE, and/or the simulations of the initial lattice model. 
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